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The next generation of weak lensing surveys will trace the evolution of matter perturbations and 
gravitational potentials from the matter dominated epoch until today. Along with constraining the 
dynamics of dark energy, they will probe the relations between matter overdensities, local curvature, 
and the Newtonian potential. We work with two functions of time and scale to account for any 
modifications of these relations in the linear regime from those in the ACDM model. We perform 
a Principal Component Analysis (PCA) to find the eigenmodes and eigenvalues of these functions 
for surveys like DES and LSST. This paper builds on and significantly extends the PCA analysis of 
Zhao et al. (2009) [1] in several ways. In particular, we consider the impact of some of the systematic 
effects expected in weak lensing surveys. We also present the PCA in terms of other choices of the 
two functions needed to parameterize modified growth on linear scales, and discuss their merits. We 
analyze the degeneracy between the modified growth functions and other cosmological parameters, 
paying special attention to the effective equation of state w(z). Finally, we demonstrate the utility 
of the PCA as an efficient data compression stage which enables one to easily derive constraints on 
parameters of specific models without recalculating Fisher matrices from scratch. 



I. INTRODUCTION 

The dynamics of linear scalar perturbations in alter- 
native theories of gravity, as well as in exotic models 
of dark energy (DE) and dark matter, can differ signif- 
icantly from that predicted by the cosmological concor- 
dance model, ACDM, even when the expansion histories 
are the same [2-6]. As in [1], we will use the term modified 
growth (MG) when referring to all these models, including 
those based on General Relativity (GR). It is expected 
that ongoing and upcoming tomographic weak lensing 
(WL) surveys, such as the Dark Energy Survey (DES) [7] 
and Large Synoptic Survey Telescope (LSST) [8], com- 
bined with the CMB and SNe data, will tightly constrain 
such modifications of growth dynamics on cosmological 
scales [1, 9-15]. 

Quite generally, a theory of gravity specifies how the 
metric perturbations relate to each other and how they 
are sourced by perturbations in the stress-energy ten- 
sor. In GR these relations are given by the anisotropy 
and Poisson equations, respectively. As in [1, 9, 16-21], 
we introduce two functions of time and scale, /i(a, k) and 
7(a, /c), to allow for general departures of these equations 
from their ACDM form. By definition, these functions 
are equal to unity in ACDM, but generally have a time- 
and/or scale-dependence in alternative models of grav- 
ity and in models with clustering DE or a significant hot 
dark matter component, such as massive neutrinos. Dif- 
ferent, but in essence equivalent, parameterizations are 
used in [13, 14, 17, 18, 22-27]. 

Such parametrizations can be used to test the valid- 
ity of ACDM in a model-independent way, which is the 
main focus of this paper. If, instead, one aims to test a 
particular theory, there is no need to use these two func- 
tions, since one can derive the exact equations from the 
action and then calculate predictions for the observables 



to constrain the parameters (typically just a few) of the 
theory. However, as we will discuss in Section VII, even 
in this case it can be advantageous to "store" informa- 
tion contained in observables into intermediate functions, 
such as our ji and 7. For instance, rather than modifying 
the standard Boltzmann codes, such as CAMB [28, 29], 
differently for each specific model alternative to ACDM, 
one can modify them once to work for arbitrary \i and 
7. Then, to evaluate observables in a given theory, it re- 
mains to work out its prediction for fi and 7. These may 
be approximate expressions, valid over a limited range of 
scales, or obtained by numerically solving a smaller set 
of equations to find fi and 7. 



In [1], a Fisher forecast and Principal Component 
Analysis (PCA) were performed to find the eigenmodes 
and eigenvalues of fi and 7 for surveys like DES and 
LSST, complemented with CMB and SNe data. The 
number of well-constrained modes approximately corre- 
sponds to the number of degrees of freedom (DoF) of 
these functions that can be measured by the surveys. 
Their scale and time dependence indicate the ranges 
where the surveys will be most likely to detect deviations 
from ACDM. The aim of this paper is to provide the de- 
tails and to expand the study of [1] in several ways. We 
present a detailed study of degeneracies between the MG 
functions and other cosmological parameters, paying spe- 
cial attention to the degeneracy with the DE equation of 
state w(z). We also present the PCA of another choice of 
MG functions that helps to demonstrate the information 
content of WL surveys. We then investigate the effect 
of some of the systematics expected in WL surveys, and 
demonstrate the utility of the PCA approach as a data 
compression stage by using it to derive constraints on a 
parameter of a specific model. 
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II. THE FORMALISM 

A. Evolution of linear perturbations 

We consider linear scalar perturbations to the flat 
Friedmann- Robert son- Walker Universe, with the choice 
of the Newtonian gauge for the metric. The line element 
reads 

ds 2 = a(r) 2 [-(l + 2^)dr 2 + (1 - 2$)dx 2 ] , (1) 

where \I/ and $ represent time- and space-dependent 
scalar perturbations of the metric and r is the confor- 
mal time. We work in Fourier space, using the same 
symbols to indicate the perturbations in space and their 
Fourier transformed counterparts, i.e \I/ = \I/(a, k) and 
<3> = 3>(a, k). The same convention is applied to the mat- 
ter density contrast 5 = Sp/p and the divergence of the 
velocity field = ik^Vj. We assume that the matter 
perturbations obey the standard conservation equations, 
which for dark matter read: 

5 + = (2) 

e + ue-kH = , (3) 

where the dot denotes differentiation with respect to con- 
formal time r, and 7~L = a~ 1 da/dr. For the sake of sim- 
plicity we ignore radiation or baryonic effects, but they 
can be easily included if relevant. 

One needs two additional equations to close the system 
for the four variables ^, 5 and 6. These are normally 
provided by a theory of gravity, which prescribes how 
the two metric potentials relate to each other, and how 
they are sourced by the matter perturbations. Since we 
aim to test potential departures from ACDM, rather than 
working with a particular gravity theory, we close the 
system of equations by introducing two general functions 
of scale and time defined via: 



kH = -47ra 2 Gp(a,k)pA (5) 

where A = S + 3H6/k 2 is the comoving matter den- 
sity perturbation. Eqs. (2)- (5) form a closed system that 
can be used to calculate the evolution of perturbations 
for any given functions p and 7; they were extensively 
discussed in [16]. There, among other things, it was ex- 
plicitly shown that they respect the superhorizon consis- 
tency conditions for adiabatic perturbations [30, 31] as 
long as (k/1-L) 2 /(p~f) when (k/U) 0. In the New- 
tonian limit the functions p and 7 are related to the 
Post-Newtonian-Parameters (PPN) of the Eddington- 
Robertson-Schiff formalism [32]. Specifically, p — >• «ppn 5 
and 7 — >• — 7ppn/<^ppn where appN and 7ppn represent 
respectively the strength of gravity and the amount of 
curvature per unit mass. 

By design, we have p = 7 = 1 in the standard cos- 
mological model ACDM. Departures of p and/or 7 from 



unity can happen if, for example, DE clusters or if it 
carries a non-negligible anisotropic stress. Alternatively, 
one could have p ^ 1 due to a significant fraction of mas- 
sive neutrinos, which free stream on small scales. Finally, 
alternative gravity models generally predict scale- and 
time-dependent p and/or 7 [6]. 

As mentioned in Introduction, the main benefit of 
using these functions is that they allow for a model- 
independent test of the growth dynamics on cosmolog- 
ical scales. Any measured deviation of either p or 7 from 
unity would signal a departure from the ACDM model. 
It should be emphasized that p and 7 do not necessar- 
ily have a simple form in specific models of MG, and 
generally depend on the choice of the initial conditions 
[33-35]. For instance, in scalar-tensor models of gravity, 
the ratio of <I> and \I/ is not a fixed function of k and a. 
Instead, it is an expression that involves the time deriva- 
tives of \I/ and This means that p and 7 correspond 
to solutions of equations of motion of a theory, rather 
than being a general prediction of a theory. Neverthe- 
less, the functions p and 7, while phenomenological in 
nature, are theoretically consistent and allow us to test 
for departures from ACDM independent of how complex 
the underlying theory of gravity is. One simply needs to 
be careful when translating the findings on p and 7 into 
results on the parameters of specific models [36], paying 
attention to the choice of the initial conditions and pos- 
sible additional simplifications, such as the quasi-static 
approximation. 

Depending on the circumstances, such as the type of 
data available or the type of theory one wants to test, it 
can be more convenient to replace either p or 7 with the 
function E(a, k) defined as 

k 2 ($ + 9) ee -87ra 2 GE(a, k)pA (6) 

The advantages of using different combinations of /i, 7 
and E are discussed at length in [16]. In this paper we 
will present and compare the results for (/i, 7) and (/i, E). 

The implementation of this formalism in CAMB, which 
uses the synchronous gauge, is detailed in [9, 37]. 



B. Principal Component Analysis 

Our goal is to determine how well p(a, k) and 7(a, k) 
can be constrained by future surveys, minimizing any as- 
sumption on the functions themselves. Therefore, rather 
than employing a specific expression for p and 7, we will 
treat them as unknown functions of both time and scale, 
and bin them on a grid in the (z, k) space (notice that we 
are now using the redshift z as the time variable). With 
m z-bins and n k-bins, we have m x n grid points to 
which we associate a value of the two functions. This is a 
2 x 2-dimensional problem, and in [1] we indicated these 
values with pij and 7^. However, in practice, we trans- 
formed the 2D problem into a ID one by mapping the 
grid into a chain, therefore transforming the matrices of 
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values into two m x n-dimensional vectors. We shall indi- 
cate the components of the vectors with and ji where 
i = 1, . . . , m x n. In addition, we also bin the DE equa- 
tion of state w(z) in redshift, creating a m-dimensional 
vector and vary the usual cosmological parameters: the 
Hubble constant h, cold dark matter density Vt c h 2 , the 
baryon density Vt^h 2 , the optical depth r, the scalar spec- 
tral index n s , and the amplitude of scalar perturbations 
A s . We assume that the bias is scale-independent on the 
linear scales considered in our analysis and introduce 
constant bias parameters, one for each photometric bin 
of the survey. 

We then use the Fisher matrix formalism to esti- 
mate the anticipated covariance of our parameters pi, 
i = 1, . . . , 2xmxn+m+6+7V5. According to the Cramer- 
Rao theorem, any unbiased estimators for the parameters 
will give a covariance matrix that is not better than the 
inverse of the Fisher matrix of the parameters. There- 
fore, after choosing our observables and experiments as 
described in next Section, we build the Fisher informa- 
tion matrix for the parameters pi. Then, we invert it to 
determine the anticipated covariance matrix 

dj = ((Pi-Pi)(Pj -Pj)) , (7) 

where pi are the assumed best fit, or "fiducial", values. 

Suppose that we want to know how well a given com- 
bination of experiments will measure \i. We marginalize 
over the other parameters, and consider the \i block of 
the covariance matrix, C^. Since the individual pixels 
of fi are highly correlated, the covariance matrix will be 
non-diagonal, and the value of \i in any particular bin 
will be practically unconstrained. The PCA is a way to 
decorrelate the parameters and find their linear combi- 
nations that are best constrained by data. Namely, we 
solve an eigenvalue problem to find a matrix that 
diagonalizes C M : 

where A^'s are the eigenvalues. Smaller values of cor- 
respond to the better constrained linear combinations of 
/i's: 

mxn 

One can think of a's as the new set of uncorrelated pa- 
rameters obtained by a rotation of /i's, with the error on 
OLi given by \f\{. In practice, one finds that only a few 
of the a's are well constrained (i. e. their A's are small), 
while most are essentially unconstrained. This is the main 
benefit of performing a PCA - it takes a function with 
many (infinite) degrees of freedom and isolates their few 
uncorrelated linear combinations that can be constrained 
by a given experiment. By construction, W T W = /, so 
Eq. (9) can be inverted as 

mxn 

Mi-Mi= E W tj<*i ' ( 10 ) 



where i labels a point on the (z, k) grid and j label the 
eigenmode. Thus, taking the continuous limit, and using 
p, = 1 as the fiducial value, we can formally rewrite this 
as 

n(z,k) = l + Y j ajW?(z,k) , (11) 

3 

which is an expansion of \i into an orthogonal basis of 
eigenvectors Wj{z,k). We can now rearrange the eigen- 
vectors into a 2D form and plot them as surfaces in the 
(z, k) space. We will refer to these surfaces as the princi- 
pal components (PC's) or eigenmodes] the shapes of the 
best constrained eigenmodes indicate the kind of features 
in ji that experiments are most likely to constrain. The 
regions in (z, k) where the best constrained eigenmodes 
peak indicate the sweet spots, i.e. the intervals in time and 
scale where a given experiment will be more sensitive to 
departures from ACDM. The number of nodes in the z 
and k directions indicate the degree of sensitivity of the 
surveys to the z- and /^-dependence of fi. The same pro- 
cedure can be repeated for the function j(a : k) (or w{a)), 
in this case isolating and diagonalizing the 7 block of the 
covariance matrix to determine the eigenvectors W 7 (z, k) 
and the corresponding eigenvalues. 

The procedure outlined above addresses the ability of 
data to constrain \i and 7 separately 1 . Namely, when 
deriving the eigenmodes and eignevalues of \i (7) we 
marginalize over 7 (/i). However, observations probe com- 
binations of the potentials and ^ that depend on both 
\i and 7. This yields a degeneracy between \i and 7 and, 
by marginalizing over one, we lose the information that 
is common to both functions. In addition to forecast- 
ing separate constraints on /i and 7, one may want to 
know how sensitive data is to any departure from the 
standard growth. Namely, we may want to answer a less 
ambitious, yet equally interesting, question of whether 
either of the two functions deviates from unity, without 
specifying which. For this purpose, we want to save the 
information common to both functions, which we previ- 
ously lost by mutual marginalization. Hence, we consider 
the combined principal components of \i and 7. We follow 
the same procedure as before, except now we diagonalize 
the block of the covariance matrix containing the pixels 
of ji and 7. The components of the matrix that diagonal- 
izes C M7 will be where now z, j = 1....2m x n; each 
eigenmode j consists now of a double series of pixels on 
the (fc, z) grid, representing combined eigenmodes in the 
two sub-spaces. 



1 We note that in this paper we do not attempt to reconstruct 
fi(a, k) and 7(0, k) from data, nor we propose PCA as a recon- 
struction tool. Instead we forecast the ability of surveys to detect 
departures of 11 and 7 from unity and use PCA to determine the 
best constrained eigenmodes. 
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C. Covariance matrix for E 

As discussed in [16], the choice of two functions 
parametrizing MG is not unique. Depending on the cir- 
cumstances, it can be preferable to replace 7 with the 
function E defined in Eq. (6), and work with the com- 
bination (/i,E). In that case, one way to proceed is to 
repeat the procedure outlined in Sec. II B for the new 
combination (/i,E). Alternatively, one can use the infor- 
mation already stored in the (/i,7) pixels to infer the 
covariance matrix for E, which is what we proceed to do. 
From Eqs. (4)- (6), we have 

£ = ^(1 + 7)- (12) 

Then, pixelizing E on the same (fc, z) grid, we can derive 
its covariance matrix in terms of the covariance matrix 
elements for \i and 7 as 

= ^^ + (1 + 70(1 + 7^" 

+ ^(1+^)^ + ^(1 + 7^], (13) 

where, for example, C^J 1 is the covariance between ji and 
fij . Analogously, one can derive the covariance of E with 



III. THE OBSERVABLES 

The ongoing and future tomographic large scale struc- 
ture surveys (like DES [7], PAN-STARR [38] and 
LSST [8]) will provide measurements of galaxy number 
counts (GC) and weak lensing (WL) on the same patch 
of sky, as well as a large number of supernovae (SN). 
This, in combination with the full sky CMB data from 
WMAP and Planck, allows us to consider all possible 
cross-correlations: CMB-WL, CMB-LSS, LSS-WL, in ad- 
dition to the three auto-correlation functions. This infor- 
mation can be used to forecast constraints on the differ- 
ences between the metric potentials and the scale-time 
variation of the effective Newton constant, parametrized 
respectively by 7 (4) and \i (5). In what follows, we pro- 
vide a brief overview of the observables and the assump- 
tions about the experiments which only have minor dif- 
ferences from the assumptions used in [9]. 



A. Angular spectra 

For any two fields, X(h) and Y(h)^ measured by an 
observer as function of the direction on the sky n, one 
can define the angular power spectrum Cf Y via 

c x Y{e) = f2 2 l±lcrPe(co S e), (14) 

£=0 71 



where C XY (0) = C XY (\&! - n 2 |) = (X(n 1 )y(n 2 )> is the 
two-point correlation function, and Pi are the Legendre 
functions. In a flat universe, Cf Y can be expressed in 
terms of the primordial curvature power spectrum 
and the angular transfer functions lf lY (k) as 

cr=^Jf^W(k)ir(k), as) 

where Ig are the transfer functions defined as 

If (k) = dzW x {z)ji[kr(z)]X{k,z). (16) 
Jo 

and similarly for I Y . In the above, z* is a sufficiently 
high redshift at which the initial condition for the mode 
k is specified, ji are the spherical Bessel functions, r(z) 
the comoving distance to a point at redshift z, and 
Wx are the window functions which, depending on the 
observable, specify the range of redshifts contributing 
to X. Finally, X(k,z) is the Fourier transform of the 
three-dimensional field X{nr(z),z) responsible for pro- 
ducing the two-dimensional observable X (i.e. X(h) = 
J °° dz Wx(z)X(nr(z), z) ). A detailed derivation of the 
above expressions is given in [9]. We adopt adiabatic 
initial conditions as detailed in [39, 40]. The observable 
quantities for which we evaluate J^'s are GC in several 
photometric redshift bins, the WL shear in several bins, 
and the CMB temperature anisotropy. 
For GC, we have 

I?*(k) = b % dzW Gz (z)j £ [kr(z)}5(k,z) , (17) 
J 

where bi is the bias, WqXz) is the normalized selection 
function for the zth redshift bin, and <S(fc, z) is the density 
contrast transfer function (and we have dropped the tilde 
following the convention outlined in Sec. II A). We work 
under the assumption that on large scales the bias can 
be treated as scale-independent and can be modeled with 
one free parameter bi for each redshift bin i. 
For weak lensing, the relevant Jg's are given by 

if* (k) = dzW K% (z)ji [kr{z)\ (* + $) , (18) 
Jo 

where W K .(z) is the window function for the ith bin of 
sheared galaxies with a normalized redshift distribution 
W Si (z), i.e.: 

wm = r dz' r{z,) ;^ z) ws^') • (19) 

The transfer functions for the CMB temperature 
anisotropy receive contributions from the last-scattering 
surface (at z ~ 1100) and from more recent redshifts via 
the Integrated Sachs- Wolfe (ISW) effect. The modifica- 
tions of gravity considered in this paper are negligible at 
recombination. Therefore, their only imprint on the CMB 
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will be via the ISW effect. For the ISW contribution to 
the CMB, we have 



tISW 



(k) 



dze- T ^j £ [kr(z)}^-[^^^} , (20) 
o az 



where r(z) is the opaqueness function. 

We numerically evaluate the transfer functions lg us- 
ing the first version of the publicly available code MG- 
CAMB (Modified Growth with CAME) [9] 2 and obtain 
the angular spectra C* Y . A joint analysis of CMB 
and data from a tomographic lensing survey with M 
GC redshift bins and N WL bins gives us a total of 
3 + M(M + l)/2 + N(N + l)/2 + M + N + MN dif- 
ferent types of CVs respectively from CMB, GC, WL, 
GCxCMB, WLxCMB and WLxGC, (we do not corre- 
late CMB polarization with GC and WL). For example, 
combining Planck with DES, with M = 4 GC bins and 
N = 4 WL bins, gives us 47 different types of spectra. A 
combination of Planck with LSST, with 10 GC bins and 
6 WL bins, gives us 155 different C/s. 

We only use the parts of the spectra that correspond 
to the linear cosmological regime. Including higher £, or 
smaller scales, would require us to account for non-linear 
effects which, strictly speaking, are not allowed within 
our framework. To accurately model growth on non- 
linear scales, one needs input from N-body simulations, 
which can only be performed for specific modified gravity 
theories. The fact that we are not testing a specific model, 
but constraining a general departure from GR, defined in 
terms of linear perturbation variables, precludes us from 
having a reliable description of non-linear corrections. We 
restrict ourselves to the linear regime by cutting off the 
Cf Y spectra at £ max ~ 0.2 hr(z s ). This cutoff roughly 
corresponds to /c max ~ 0.2/iMpc -1 at z = 0. There is 
certainly a wealth of information about MG parameters 
on smaller scales [41], and while it would be tempting to 
include it in our analysis, it would make our predictions 
unreliable since our analysis is limited to linear theory. 



B. Fisher matrices 

To determine how well the surveys will be able to con- 
strain our model parameters, we employ the standard 
Fisher matrix technique [42]. The inverse of the Fisher 
matrix F a b provides a lower bound on the covariance ma- 
trix of the model parameters via the Cramer- Rao inequal- 
ity, C > F _1 [42]. For zero-mean Gaussian-distributed 



2 A new version of MGCAMB was recently introduced in [37], 
and is publicly available at http : //www . sf u . c a/ ~aha25 /MGCAMB . 
html. The first version was based on CAMB-Sources, which made 
it easy to evaluate WL and GC spectra, but was not compati- 
ble with CosmoMC. The later version is written as a patch for 
CAMB, is compatible with CosmoMC, but does not evaluate WL 
and GC spectra yet. 



observables, like the angular correlations Cf Y , the Fisher 
matrix is given by 



F a 



/sky Yl 

£=£mi: 



21- 



dp a dp b 



, (21) 



where p a is the a th parameter of our model and is the 
"observed" covariance matrix with elements Cf Y that 
include contributions from noise: 



O p — w p 



N, 



XY 



(22) 



Eq. (21) assumes that all fields X(n) are measured over 
contiguous regions covering a fraction f s ^ y of the sky. 
The value of the lowest multipole can be approximately 
inferred from £ min « n/ (2/ s k y ). It is also possible to write 
expressions for separate contributions to the Fisher ma- 
trix from particular subsets of observables, as detailed in 
[9]. 

The noise matrix Nf Y includes the statistical noise as 
well as the expected systematic errors. Systematics are 
notoriously difficult to predict, and are often ignored in 
forecasts. In this paper, we consider two cases - the case 
when N* Y includes the statistical noise only, and the 
case when certain types of systematics are included. We 
follow [43, 44] and consider three sources of systematics 
for future tomographic surveys: photo- z errors, as well 
as additive and multiplicative errors due to the uncer- 
tainty of the point spread function (PSF) measurements. 
The details of our modeling of the systematic effects are 
presented in Appendix A. Our assumptions about the 
statistical errors in CMB, CG and WL are presented in 
the following subsection. 

For supernovae, the information matrix is 



• ab 



N 



1 dm(zi) dm(zi) 
o-(zi) 2 dp a dp b 



(23) 



where m(z) is their redshift-dependent magnitude, the 
summation is over the redshift bins, and a(zi) = 0.13 
(see Sec. IIIC2 for details). 

Given a set of theoretical covariance matrices over a 
given multipole range, and the specifications for the ex- 
pected noise in specific experiments, we can compute the 
Fisher matrix. The derivatives with respect to the pa- 
rameters (pixels of \i and 7 in our case) are computed 
numerically using finite differences. 



C. Experiments 

The data considered in our forecasts include CMB 
temperature and polarization (T and E), WL of distant 
galaxies, GC, their cross-correlations, and SNe observa- 
tions. We assume CMB T and E data from the Planck 
satellite [47] , the galaxy catalogues and WL data by the 
Dark Energy Survey [7] and Large Synoptic Survey Tele- 
scope (LSST) [8], complemented by a SNe data set pro- 
vided by the futuristic Euclid-like survey [48]. 
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1. DES and LSST 



We take the total galaxy number density to be given 



N G (z) oc z 2 exp(-(z/z ) 2 ) , 



(24) 



which is a slight modification of the model due to 
Wittman et al. [49]. The parameter zq depends on the 
experiment and defines the redshift at which the most 
galaxies will be observed. The galaxies can be divided 
into photometric redshift bins, labeled with index i, 



N G {z) = Y^N G ^)- 



(25) 



In our analysis, we assume that the photometric redshift 
errors are Gaussian distributed, and that their rms fluc- 
tuations increase with redshift as a(z) = cr max (l+2:)/(l + 
z max ). The bin sizes are chosen to increase proportionally 
to the errors. The resulting photometric redshift distri- 
butions are given by 



1 



N Gi (z) = -N G (z) 



erfc ( 



Zi-l 



V V2cr(z) 



erfc ( 



(26) 

where erfc is the complementary error function. For a 
given photometric redshift bin, the normalized window 
function that appears in Eq. (17) is given by 



W Gi (z) 



N Gl (z) 



(27) 



where N % is the total number of galaxies in the ith. bin. 

DES is a project aimed at studying the nature of 
the cosmic acceleration, and is planned to start obser- 
vations in 2012 [7]. DES includes a 5000 square degree 
multi-band, optical survey probing the redshift range 
0.1 < z < 1.3 with a median redshift of z = 0.7 and 
an approximate l-cr error of 0.05 in photometric red- 
shift. In our simulation, for both WL and GC, we as- 
sume a sky fraction / s k y = 0.13, and an angular den- 
sity of galaxies N G = 10 gal/arcmin 2 . We also assume 
7rms = 0-18 + 0.042 z, which is the rms shear stemming 
from the intrinsic ellipticity of the galaxies and mea- 
surement noise, and the photometric redshift uncertainty 
given by a(z) = 0.05(1 + z). 

LSST is a proposed large aperture, ground-based, wide 
field survey telescope [8]. It is expected to cover up to 
half of the sky and catalogue several billion galaxies out 
to redshift z ~ 3. For LSST forecasts, we adopt param- 
eters from the recent review paper by the LSST col- 
laboration [50]. Namely, we use / s k y = 0.5, N G = 50 
gal/arcmin 2 for both WL and counts, 7 rms = 0.18 + 
0.042 z, and a(z) = 0.03 (1 + z). 

For both DES and LSST, we take the GC photometric 
bins to be separated by 5a(z). This leads to four redshift 
bins for DES and ten for LSST. For WL (source) galaxies, 
we use four bins for DES and six for LSST. 



2. CMB, Supernovae, and other priors on cosmological 
parameters 

In our forecasts, we assume spatially flat geometry and 
parametrize the dark energy equation of state using 20 
bins from z = to z = 3 uniform in z, and one wide bin 
from z = 3 to z = 1100. In addition to the MG parame- 
ters, and the = M bias parameters, we vary h, Vt c h? , 
Vtbh 2 , r, n s , and A s . Their fiducial values are taken to be 
those from the WMAP 7-year data best fit [51]: Sl h h 2 = 
0.023, n c h 2 = 0.11, h = 0.71, r = 0.088, n s = 0.963. The 
fiducial values for bias parameters are motivated by the 
parametrized halo model described in [52]. Imposing a 
prior on the value of h from the Hubble Space Telescope 
(HST) did not make a noticeable difference in our results. 

We assume the expected CMB data from the Planck 
mission [47] of the European Space Agency (ESA) us- 
ing the same parameters as in [9]. In addition, to bet- 
ter constrain the background expansion parameters, we 
include simulated SNe luminosity data for a Euclid-like 
survey, e.g., the one proposed in [48]. We generate 4012 
data points randomly distributed in 14 redshift bins from 
z = 0.15 to z = 1.55, and combine 300 low-z SNe from 
the Nearby Supernova Factory (NSNF) survey [53]. We 
calculate the exact distance modulus for each model, and 
put a Gaussian noise with rms error of a = 0.13 to dis- 
place all the data points. The absolute magnitude, or the 
so-called nuisance parameter M., is treated as an unde- 
termined parameter in our analysis. 

Note that DES and LSST will produce their own SNe 
luminosity measurements that we did not include in our 
forecasts. Instead, we assume that a high quality SNe 
catalogue will eventually become available, and focus on 
the dependence of MG constraints on the quality of the 
WL and GC data. 



IV. THE PCA OF LINEARIZED MG 

In this Section we analyze the principal components 
(eigenmodes) of the functions \i and 7 defined in Sec. II 
for the combination of large scale structure (WL and 
GC), CMB and SNe experiments detailed in Sec. IIIC. 
In particular, we want to investigate the effect of degen- 
eracies with other cosmological parameters on the eigen- 
modes and eigenvalues of \i and 7. For this purpose, in 
Subsection IV A, we start with the simplest case where 
only uncertainties in the pixelated functions \i and 7 (re- 
ferred to as "MG pixels") are considered, with all other 
parameters fixed to their fiducial values. The highest red- 
shift pixels (3 < z < 30) of /i and 7 are outside the range 
directly probed by the WL surveys. However, they do 
impact the observables: /i(z > 3, k) re-sets the amplitude 
of the growth at all lower redshifts, while j(z > 3, k) 
affects the ISW contribution to the CMB temperature 
anisotropy. The sensitivity of observables to variations in 
/i(z > 3, k) and 7(2 > 3, k) depends on the assumed high- 
z cutoff - making the high-z interval wider increases the 
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FIG. 1: The eigenmodes of \i and 7 for LSST(+SN+CMB) 
with all other parameters fixed to fiducial values. 



FIG. 2: The combined eigenmodes of (m, 7) for 
LSST(+SN+CMB) with all other parameters fixed to 
fiducial values. 



sensitivity. Thus below we focus on quantities that are 
independent of this cutoff, such as the eigenmodes and 
eigenvalues of \i and 7 obtained after marginalizing over 
the high-z MG pixels, as well as the uncertainty in the 
overall growth between recombination and z = 3, which 
is directly controlled by /i(z > 3, fc), but is independent 
of the choice of the high-z cutoff (see IV D). In IV B, 
we discuss effects of marginalizing over the high-z pix- 
els and galaxy bias parameters. In (IV C), we consider 
degeneracies with the cosmological vanilla parameters 
{Qbh 2 , Q c h 2 , h, r, n s , A s } and the effective equation of 
state w(z). Throughout this Section, we always marginal- 
ize over the SNe nuisance parameter. To gain additional 
insight, we also project our findings onto the function E. 

We examine the parameter degeneracies using LSST 
as the fiducial survey for WL and GC. Then, in IV E, we 
compare the final case, in which we marginalize over all 
the parameters except the \i and 7 in the range < z < 3, 
to the analogous forecast for DES. 



A. PCA for LSST with all parameters, except MG 
pixels, fixed (high-z information included) 

As a first step, we study the eigenmodes and associated 
uncertainties of the MG functions without considering 
their covariance with other parameters. We do this by 
isolating and inverting the block of the Fisher matrix 
containing only MG pixels (including the high-z bins) 
and diagonalizing the resulting covariance matrix. 

Fig. 1 shows some of the eigenmodes of \i and 7. Each 
panel in these plots represents a region in (z, k) space 
with an eigenmode function plotted as a surface, as de- 
scribed in Sec. II A. We order, and consequently number, 
the eigenmodes according to how well they can be con- 
strained, i.e. following the ordering of the corresponding 
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FIG. 3: The eigenmodes of E for LSST(+SN+CMB) with all 
other parameters fixed to fiducial values. 
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FIG. 4: The uncertainties (square roots of covariance eigneval- 
ues) associated with the eigenmodes of /x, 7, E and the com- 
bined (/i, 7) case for LSST(+SN+CMB) with all other param- 
eters fixed to fiducial values. 



8 



errors (square root of the covariance matrix eigenvalues) 
from the smallest (best constrained) to the largest (least 
constrained) . 

The first feature to notice is in the /i eigenmodes (top 
two rows in Fig. 1), where the best constrained modes 
peak at high-z, and show no features at low z; we need 
to get to the 16th mode to start seeing some features in 
the z < 3 interval, which is the actual redshift range of 
LSST. This is because fi directly affects the growth of 
matter density perturbations, and changing the ampli- 
tude of perturbations at some redshift affects the growth 
at lower redshift s. As we discuss in IV D, the constraint 
on fi at z > 3 depends strongly on the width of the bin, 
i.e. the value of the arbitrarily chosen high- z cutoff. It is 
also highly correlated with \i in lower redshift bins and 
some of the vanilla cosmological parameters. 

Unlike /i, the best constrained modes of 7 do not have 
a support at high z. They peak at low redshifts (bottom 
two rows in Fig. 1). This follows from the fact that, ac- 
cording to our definition, to measure 7 one needs to mea- 
sure both and Therefore, bounds on 7 come mostly 
from combining the GC data, which probes ^ (affected 
by /i), with the WL signal probing (<£ + \I>) (affected by 
both \i and 7). One can also see that there are modes 
that peak at z > 3 and at low k (very large scales). This 
is mainly due to the ISW effect, seen as a contribution 
to the CMB temperature spectrum, which is sensitive to 
the time variation of both potentials at all times after 
the last-scattering. 

An important thing to notice is that the best con- 
strained modes of \i and 7 show oscillations (nodes) in 
/c, but no nodes in the z-direction. One has to look at 
the higher number modes to start seeing oscillations in 
z. For instance, for /i, the first z-node appears in the 
26th eigenmode. The number of nodes is indicative of 
the sensitivity of the function to changes in k and and 
we see that the experiments are significantly more sen- 
sitive to scale-dependent features of the MG functions, 
and not as sensitive to the time-dependence. As already 
pointed out in [1], this is expected, since the impact of a 
scale-dependent change in fi (or 7) is directly translated 
into a scale-dependent feature in the WL and GC spec- 
tra. For instance, in the case of GC, \i effectively appears 
as a scale-dependent normalization factor. On the other 
hand, the projection of time-dependent features of the 
MG functions onto the observables involves integration 
over time which makes detecting z-dependent features 
harder. Also, the amount of information coming from the 
radial (z-direction) is limited by the fact that LSST only 
probes structures at z < 3, and by the fact that we con- 
sider only linear scales, effectively cutting off a significant 
volume at low z. 

Fig. 2 shows the best combined eigenmodes of (/i, 7). 
Every combined eigenmode is represented by a pair of 
values at each point on the (fc, z) grid, one value result- 
ing from a variation of \i on that grid point and the other 
from a variation of 7 on the same grid point. We show 
each eigenmode as a pair of surfaces, one corresponding 



to ji and the other to 7. It can be noticed that the best 
constrained modes peak at high redshift, which is the re- 
sult of high sensitivity of the growth to changes in \i in the 
high-z bins. Since 7 does not directly affect the growth 
rate, the 7 surfaces of these eigenmodes have a very low 
amplitude, requiring us to amplify them by a factor of 100 
in order to make them visible in the plots. The combined 
eigenmodes that peak at low redshifts, starting from the 
13th, do not exhibit separation of scale or time dependent 
oscillations, but rather have a diagonal form in the (z, k) 
space, showing a degeneracy between scale and time. This 
is because the WL observables dominate the information 
for combined modes at low redshift. Indeed, the changes 
in the weak lensing kernel due to a shift of the lens along 
the line of sight (i.e. a change in redshift) are degenerate 
with those due to a resizing of the lens (i.e. a change in 
scale) . 

Fig. 3 shows the best eigenmodes for the function E. 
As it is clear from its definition (6), E is directly sensi- 
tive to the lensing potential ($ + therefore its signal 
comes mainly from low redshifts and it does not have 
well constrained high redshift modes. For this reason the 
plots for E span only over the low redshift interval. From 
these plots we notice a k — z degeneracy analogous to 
the one found in the combined (/i, 7) modes as it is clear 
from a comparison of Fig. 2 and Fig. 3. 

Finally, Fig. 4 shows the uncertainties (square roots of 
covariance eigenvalues) associated with the eigenmodes 
of /i, 7, E and the combined modes. As expected, the 
combined eigenmodes are constrained best, since they 
contain all of the information about any departure of 
either \i or 7 from unity. The constraints on 7 are the 
weakest, since 7 is not directly constrained by WL nor 
GC, while \i does better than 7 since it is directly con- 
strained by GC. The eigenmodes of \i that peak at z > 3 
are constrained better than those of E, while E is mea- 
sured better than ji for modes with support at lower red- 
shifts. This is because at low redshifts most of the in- 
formation comes from WL, which is directly sensitive to 
E. As mentioned earlier, the bounds on the z > 3 bin of 
fi are dependent on the arbitrary cutoff and in the next 
subsection we will marginalize over it. 



B. PCA for LSST after marginalizing over high z 
and galaxy bias, everything else is fixed 

In order to remove the dependence of our results on 
the arbitrary choice of the upper cutoff of the high-z 
bins (3 < z < 30), we marginalize over the high-z pix- 
els. This is achieved by removing the rows and columns 
corresponding to high-z pixels from the covariance ma- 
trix of the previous subsection, and then diagonalizing it 
to find the eigenmodes and eigenvalues. This essentially 
removes all the information about the growth at z > 3 
and, with that, all the z > 3 features in the eigenmodes. 
This is seen in Figs. 5-7, where the eigenmodes for /i, 7, 
the combined (/i, 7) and E are shown for this case. 
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FIG. 5: The eigenmodes of fi and 7 for LSST(+SN+CMB) 
after marginalizing over the high-z bins, with all other pa- 
rameters fixed to fiducial values. 
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FIG. 6: The combined eigenmodes of (m, 7) for 
LSST(+SN+CMB) after marginalizing over the high-z 
bins, with all other parameters fixed to fiducial values. 



From Fig. 5 we notice that, in the absence of the high- 
z information, the first few best constrained eigenmodes 
of \i and 7 have similar shapes. However, this similar- 
ity fades for the higher order eigenmodes, reflecting the 
different dependences of the two functions on the metric 
potentials. For instance, the first node in z appears at 
the 9th mode for \i and only at the 12th for 7, reflecting 
a higher sensitivity of \i to time dependent features. 

The best constrained combined (/i, 7) eigenmodes have 
the same shapes in the /> and 7-surfaces, but the 7- 
surfaces have a lower amplitude (Fig. 6). This is again 
explained by the fact that both functions are constrained 
by the same experiments (WL and GC), therefore having 
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FIG. 7: The eigenmodes of S for LSST(+SN+CMB) after 
marginalizing over the high-z bins, with all other parameters 
fixed to fiducial values. 



similar eigenmodes, with \i being more directly related to 
the growth of structure. 

As can be observed from Fig. 7, after marginalizing 
over the high-z bins, the best eigenmodes of E are very 
similar to the combined (/i, 7) modes. Indeed, once the 
high-z information is removed, the combined modes are 
primarily constrained by WL and, therefore, carry more 
or less the same information as the E modes. 
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FIG. 8: The uncertainties associated with the eigen- 
modes of /i, 7, E and the combined (/i, 7) modes for 
LSST(+SN+CMB), after marginalizing over the high-z bins 
only (M.o. high-z), and after marginalization over the high-z 
bins and the galaxy bias parameters (M.o. high-z+bias), with 
all other parameters fixed to fiducial values. 

Fig. 8 shows the uncertainties associated with the 
eigenmodes of /i, 7, E and the combined (/i,7) modes 
after marginalizing over the high-z bins. Comparing to 
Fig. 4, we see that now there is no crossing over of the 
errors on \i and E, i.e. all the PCs of E do better than 
the fi ones. This is due to the disappearance of the eigen- 
modes with high-z support. In addition, we notice a small 
overall degradation of constraints, which is due to throw- 
ing away the information common with the low-z bins. 
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FIG. 9: The eigenmodes of fi and 7 for LSST(+SN+CMB), 
after marginalizing over the high-z bins and the galaxy bias 
parameters, with all other parameters fixed to fiducial values. 



Next, we marginalize over the galaxy bias parameters. 
This is achieved by inverting the part of the Fisher ma- 
trix that includes the MG pixels and the bias parame- 
ters. This effectively removes the information about the 
overall normalization of fi. This is manifested in the dis- 
appearance of some of the eigenmodes, including the ho- 
mogeneous eigenmodes, e.g. present in Fig. 5 but not in 
Fig. 9. The same happens to the eigenmodes of £ and of 
the combination (/i, 7) and we do not plot them. 



C. PCA for LSST after marginalizing over the 
vanilla cosmological parameters and w(z). 

We now marginalize over the cosmological parameters: 
{Q c h 2 , Qbh 2 , h, r, n s , A s } and the binned values of w(z). 
By doing so we account for the covariance of MG pix- 
els with the vanilla ACDM parameters and the effective 
dark energy equation of state. Figs. 10-13 show the PCA 
results for this case. 

The impact of marginalizing over w(z), as opposed to 
setting it to w = —1, is not dramatic and we do not 
separately show the eigenmodes for the latter case. The 
associated eignevalues, plotted in Fig. 13, show only a mi- 
nor differences. This is, in part, due to the high quality 
of the assumed SNe dataset (see Sec. IIIC 2). However, it 
is also because of marginalizing over the vanilla param- 
eters, the galaxy bias parameters and the high-z bins, 
which already throws away most of the information that 
is common between the MG pixels and w(z). The effect 
of marginalizing over bias parameters and w bins is mi- 
nor changes in the shape of modes (e.g. second mode 
in Fig. 10 to be compared with second mode in Fig. 9), 
disappearance of some of the modes and an overall degra- 
dation of constraints. 



FIG. 10: The eigenmodes of fi and 7 for LSST(+SN+CMB) 
after marginalizing over all other parameters. 
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FIG. 11: The combined eigenmodes of (m, 7) for 

LSST(+SN+CMB) after marginalizing over all other 
parameters. 
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FIG. 13: The uncertainties associated with the eigenmodes of 
/x, 7, E and the combined (^,7) modes for LSST+SN+CMB 
after marginalizing over all other parameters. Two cases are 
shown: with w bins fixed to their fiducial value of —1 (dots 
joined by lines), and with w bins varied and marginalized over 
(triangles) . 



D. Constraints on the growth at high redshift 



r (C rr = y/<rf) from 

Edr 
g-°Pi , ( 29 ) 

i 

where p^s are all the parameters of our model that affect 
r. Our calculation shows that <j r /r > 1 for all k. In other 
words, r is completely unbounded, as expected. Note that 
the calculation would need to be adjusted for DES, for 
which the highest redshift is z = 1.3, instead of z = 3. 

E. Comparison with DES 

To get a further insight into how the MG con- 
straints depend on the properties of the experiments, 
we compare PCA for LSST+CMB+SNe with that for 
DES+CMB+SN. That is, we interchange LSST with 
DES for the WL and GC data, while keeping the assump- 
tions for the CMB and SNe data the same. We perform 
this comparison only for the case in which we marginal- 
ize over all parameters, i.e. high-z bins, bias parameters, 
cosmological vanilla parameters and w(z), as in IV C. 

Figs. 14 and 15 show the eigenmodes of /i, 7 and (/i, 7) 
for DES+CMB+SN. These results should be compared 
to Figs. 10 and 11 for LSST. We choose to show a smaller 
number of eigenmodes, since there are not as many well 
constrained modes for DES as there are for LSST. 



The growth at z > 3 is not directly probed by the large 
scale structure surveys we are considering in this paper. 
Any modification to growth at z > 3, such as due to 
variations of high-z MG pixels, is observed as an overall 
shift of the amplitude of the growth at all subsequent (i.e. 
lower) redshifts. This can, in principle, be compensated 
by counter- variations of the low- z MG pixels. However, 
we find that bounds on high-z pixels are still relatively 
tight (at a percent level), even after marginalizing over 
the low-z bins. Such tight constraints are due to the large 
width of high-z pixels. Namely, a small change in value 
of \i in the high-z bins results in an accumulated modifi- 
cation of growth that can only be compensated by a very 
large variation of the low-z pixels. On the other hand, 
the low-z pixels are directly constrained by the surveys 
and large variations are not allowed. Of course, one can 
always make the bounds on the high-z bins arbitrarily 
weak by making the bins narrower. 

Since the width of high-z pixels is a somewhat arbi- 
trary parameter, one can ask if another quantity can be 
introduced to quantify the growth at z = 3. We take this 
quantity to be the ratio r(k), defined as: 



DES+SNe+CMB 



ik) 



A(z = 3,k) 

A(z rec , k) 



(28) 



where A is defined in Eq. (5) and z rec denotes the red- 
shift at recombination. We can estimate the variance in 
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FIG. 14: The eigenmodes of /x and 7 for DES+SN+CMB after 
marginalizing over all other parameters. 

As can be seen from Figs. 14 and 15, the range over 
which the eigenmodes vary is limited to smaller redshifts 
(z < 1), which reflects the redshift range probed by DES. 
Like in the case of LSST, there is a higher sensitivity 
to scale-dependent features. Furthermore, in the case of 
DES, there are modes with oscillations in z. Some of the 
combined eigenmodes are also absent for DES . Fig. 16 
shows a comparison of the uncertainties associated with 
the modes for LSST and DES. 

While the overall sensitivity of DES is less than LSST 
by a factor of few, DES, when combined with CMB and 
SNe data, is still able to constrain several eigenmodes 
with better than 10% accuracy. We may not be able to 
detect a time-dependent MG feature with high confidence 
from DES, but it is possible to constrain scale-dependent 
features. 
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FIG. 15: The eigenmodes of the combined (/i,7) modes for 
DES+SN+CMB after marginalizing over all other parame- 
ters. 
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FIG. 16: The uncertainties associated with the eigemodes of 
/x, 7, E and combined (^,7) modes for DES(+SN+CMB) and 
for LSST(+SN+CMB) errors, after marginalization over all 
other parameters. 



apparent, but not disastrous. It is true that the catas- 
trophic photo- z errors (CPE) for WL surveys [54] can 
lead to disastrous degradation on the constraints of the 
cosmological parameters, but as shown in [54], an addi- 
tional 30,000 spectroscopic redshifts can help to control 
the bias in cosmological parameters due to CPE under 
the level of statistical errors. The most significant effect is 
to reduce the ability of LSST to detect the z-dependence 
of 11. This is a preliminary analysis, which does not in- 
clude scale-dependent systematics, which can be partic- 
ularly important in MG studies (as was found in [11] 
where PCA was applied to a set of existing data includ- 
ing CFHTLS). 
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FIG. 17: Top: The uncertainties associated with the combined 
(/i, 7) eigenmodes for LSST (left) and DES (right) for four 
cases: (I) without systematics and w fixed, (II) without sys- 
tematics but w varied, (III) with systematics but w fixed and 
(IV) with systematics and w varied. Bottom: The uncertain- 
ties for LSST (left) and DES (right) relative to Case 1. 



V. EFFECTS OF SYSTEMATIC ERRORS 

Systematics are notoriously hard to predict, and fore- 
cast results that include systematics are in general sen- 
sitive to the modeling and the assumed priors [43-46]. 
Here we consider some of the systematics that will af- 
fect future tomographic surveys and study their impact 
on our PCA results. This can give us a general insight 
on how systematic errors could affect our inferences of 
the eigenmodes and the corresponding uncertainties. The 
systematics we consider here are the photo- z errors and 
some of the errors in the measurement of the point spread 
function (PSF). These errors are modeled in [43, 44] and 
we use their parametrization for our Fisher analysis. 

The effects of the systematics are detailed in the Ap- 
pendix, and with these assumptions we study the degra- 
dation of ability of DES and LSST to constrain MG by 
marginalizing over the systematics parameters without 
applying any prior. As we show later, the degradation is 



Let us look at the errors for combined eigenmodes 
shown in Fig. 17 for the following cases: (I) no systemat- 
ics, with w fixed; (II) no systematics, with w varied; (III) 
with systematics, and w fixed; (IV) with systematics, and 
w varied. We can see that allowing for systematics de- 
grades the constraints more than allowing for variations 
in w. Also, once the systematic errors are included, al- 
lowing w to vary does not degrade the constraint further. 
This means that the MG pixels are more degenerate with 
the systematics parameters than the w bins, implying 
that the uncertainty in the galaxy distributions, which 
is basically photo-z error from the systematics, can af- 
fect the constraint on the growth more significantly than 
w does. In general, the errors on the best constrained 
modes are degraded by < 10% for both LSST and DES. 
Another interesting observation is that the largest degra- 
dation does not happen for the first few modes, but for 
the intermediate modes. This is reasonable - the first 
few modes do not have nodes in z and thus are relatively 
immune to the systematics dominated by the photo-z er- 
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FIG. 18: Top: The uncertainties associated with eigenmodes 
of ii for LSST (left) and DES (right) for four cases: (I) with- 
out systematics and w fixed, (II) without systematics but w 
varied, (III) with systematics but w fixed and (IV) with sys- 
tematics and w varied. Bottom: The uncertainties for LSST 
(left) and DES (right) relative to Case 1. 



rors. 

We have shown the errors on the eigenmodes for \i in 
Fig. 18 where the same cases as above are considered. One 
notices that degradation on /i is less than what we had 
for combined case, especially no significant changes for 
DES are found. For 7, degradation is very small simply 
because the constraints on 7 eigenmodes are very weak 
in the first place. We are not therefore showing 7 errors 
here. 

Another observation is that systematics can create 
new, or destroy existing modes, so that the modes with 
the same order in the PCA sequence in the cases with 
and without systematics can be different modes. One 
can see this clearly by looking at the eigensurfaces in 
Fig. 19, where we illustrate the three \i modes and one 
7 mode with and without systematics. They look sim- 
ilar except that the modes with systematics in general 
have more nodes in /c, indicating that systematics do not 
just dilute the constraints on the old modes, but also 
make some modes disappear (or make them very poorly 
constrained). The general trend is that eigenmodes with 
very high frequency features in k are no longer well- 
constrained after inclusion of systematics. 

Overall, we find that the inclusion of systematics re- 
sults in a noticeable, but not dramatic, dilution of the 
constraints on MG from DES. This is because photo- z 
errors would most immediately affect the z-dependence 
of MG, to which DES was already mostly insensitive even 
without systematics. The main constraints from DES, as 
can be seen from the shape of the eigenmodes, will be on 
the scale-dependence of /i and 7, and that information is 
somewhat reduced, but mostly preserved. The impact of 
the systematics on LSST could be more significant, sim- 
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FIG. 19: The eigenmodes for LSST(+SN+CMB). The up- 
per (lower) panel corresponds to the case without (with) sys- 
tematics. The 10th, 11th and 12th modes of \i without system- 
atics, are compared to the 8th, 11th and 12th modes of \i with 
systematics respectively. These modes are chosen since they 
correspond to the last z-independent and first z-dependent 
mode respectively, in the two case (with and without system- 
atics). Analogously, the 4th eigenmode of 7 without system- 
atics is compared to the 4th mode of 7 with systematics. It 
is an illustration of how LSST eigenmodes are distorted as a 
result of accounting for systematics. 
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with systematics. It is an illustration of how DES eigenmodes 
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ply because LSST has a higher potential for resolving z- 
dependent features. In this preliminary analysis, we find 
that allowing for systematic errors under the assumptions 
of [43, 44] preserves most of the scale-dependent infor- 
mation from LSST, but decreases our ability to measure 
eigenmodes of /i with z-dependent features. 



VI. DEGENERACY BETWEEN MG 
PARAMETERS AND DARK ENERGY EOS 

In Sec. IV C, we marginalized over the binned equation 
of state and analyzed the impact of this marginalization 
on the eigenmodes of /i and 7. Here we investigate how 
marginalizing over MG pixels affects the eigenmodes and 
eigenvalues of w(z). In this case, we diagonalize the block 
of the covariance matrix containing the 21 lu-bins. Fig. 21 
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FIG. 21: Top: First three eigenmodes of w(z) with (left) and 
without (right) MG included. Bottom: Errors for w(z) eigen- 
modes with and without MG included (there is a prior of 0.5 
on w bins). 



(top) shows the best three eigenmodes of w in the case 
when fx and 7 pixels are co-varied, and when they are 
fixed to their fiducial values. 

The most noticeable effect of the marginalization is a 
reduction of the amplitude of the best constrained modes 
at higher z. In other words, letting MG parameters vary 
squeezes the best constrained eigenmodes of w(z) to- 
wards low redshift. This is expected, since most of the 
information on w(z) comes from the SNe and from the 
probes of LSS at higher redshift. The latter is largely 
degenerate with MG, and so the high-z information on 
w(z) is erased. On the other hand, LSS does not con- 
tribute much to the low-z information, since we restrict 
to modes in the linear regime, therefore the marginal- 
ization has little impact on the low-redshift parts of the 
w(z) modes. 

Overall the effects of the marginalization over the MG 
functions are not dramatic and future surveys will have 
the ability to measure both w and MG. Fig. 21 (bottom) 
shows the degradation of errors on the w eigenmodes af- 
ter marginalizing over other parameters, including MG. 
Note that a prior of <Ji = 0.5 has been put on the bins of 
w . 



VII. PROJECTED ERRORS ON PARAMETERS 
OF SPECIFIC MODELS FROM PRINCIPAL 
COMPONENTS 

MG pixels can be viewed as a compressed form for 
the information we get from observation about the linear 
growth. One can perform a PC A and store this informa- 
tion in the eigenmodes and their eigenvalues. The advan- 
tage of PCA is that the information can be compressed 
and used more efficiently, namely, the well-constrained 



eigenmodes usually carry almost all of the existing infor- 
mation. 

Exploiting information stored in the eigenmodes, we 
can easily emulate any other parameterization to forecast 
parameter errors without regenerating the Fisher matri- 
ces from scratch [55]. In other words, we can treat MG 
pixels as our observables and use them to calculate the 
Fisher matrices for specific model parameters. 

The Fisher matrix can be written as 



ab 



a(3 



dp a 



y a(3 



dp b 



(30) 



where O a are cosmological observables and p a are pa- 
rameters of a specific model. This can be rewritten as 



F a b 



Q p a 



^ dO a c _ 1 dO (3 
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(31) 



where Vs are MG pixels and F^ is the ij element of 
their Fisher matrix. All we need now is to expand the 
derivatives of the MG pixels with respect to a given pa- 
rameter in the eigenmode basis. That is, for each of the 



new parameters, we find the coefficients Cj 1 such that 



dp a 



(32) 



where the sum is over all the eigenmodes and e\ is value 
of the Ith. eigenmode at the iih pixel of the 2D (fc, z) grid. 
Substituting Eq. (32) in Eq. (31), we can then find the 
Fisher matrix for the parameters of our model by simple 
projection as 



F a b - C b i \ 1 j 



(33) 



where again the sum is over all the eigenmodes, and A^'s 
are the corresponding eigenvalues of the covariance ma- 
trix for MG pixels. 

We now illustrate the details of this method by apply- 
ing it to a one-parameter model which gives a good ap- 
proximation of f(R) theories in the quasi-static limit [12] 
(and which is a customized form of the more general 
parametrization introduced in [18]) 



/i(a, k) 
7(a,fc) 



1 



|A 2 k 2 a A 



1.4-10- 8 A 2 a 3 l + A 2 /c 2 a 4 
|A 2 k 2 a 4 



1 + fA 2 Pa 4 ' 



(34) 



3 Note that since we are working on a discrete (k, z) grid, our 
equations are in a discrete form. Analogous expressions could 
be written for the continuous case. For example, Eq. (32) for a 



continuous /i function would be = J^z ^ 1 e f ( a > 
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where A 2 = B c 2 /(2Hq) and c is the speed of light. The 
parameter A is the mass of the f(R) scalar degree of 
freedom today. In [10, 56], a bound of B$ < 10 -3 at 
95% confidence level was found based on the current clus- 
ter abundance data, which extends to mildly non-linear 
scales. 

Here, we forecast the constraints on Bo from LSST, 
based only on linear scales, in two ways: using a direct 
Fisher matrix calculation and a Fisher matrix projection 
described above. We choose Bq = 0, corresponding to 
GR, as the fiducial model. We use the same combination 
of future data for the direct Fisher matrix calculations as 
we did for the PC A of MG, and use the eigenmodes and 
eigenvalues of Section IV C for the Fisher matrix projec- 
tion. The derivatives on the left hand side of Eq. (32) 
can be calculated analytically from Eq. (34). Since we 
are working on a 2D grid of k and a, we use the aver- 
aged value of these derivative expressions over each pixel 
and compute the expansion coefficients of eq. (32) nu- 
merically. Table I shows the forecasted constraints on Bq 
from LSST and DES combination data obtained using 
the two methods. The results show a reasonable agree- 
ment. 



DES 


LSST 


Direct Projection 


Direct Projection 


1.5 -10" 6 2.5 -10" 6 


3.1 • 10" 7 2.4- 10" 7 



TABLE I: Error forecasts for the Bo parameter for DES and 
LSST in combination with CMB and SNe data. The results 
for Fisher matrix projection formalism used here is compared 
to the direct Fisher matrix calculation. 

One should be careful when working with projection 
method. For example, the priors used for calculating the 
covariance matrix for MG pixels (Cij in eq. (31)) should 
be the same as the priors that would be used for a direct 
fisher analysis for a model. Ideally, the fiducial models 
should also be the same in both approaches. 

As mentioned before, one of the advantages of the PC A 
approach is that we are able to compress information 
using only the best principal components. We find that in 
order to reproduce the errors shown in the "Projection" 
columns of Table I at about 95% precision we only need ~ 
25% of the eigenmodes. Such compression of information 
can be useful, given the increasing volume of cosmological 
data. 



VIII. SUMMARY AND OUTLOOK 

In this paper we have extended the principal compo- 
nent analysis first performed in [1]. As shown in previous 
sections, upcoming and future weak lensing surveys will 
provide high precision data on the relationships between 
matter over densities, curvature of space and the Newto- 
nian potential, offering an unprecedented opportunity to 
test GR on cosmological scales. 



In Section II we have introduced the MG functions 
needed to parametrize the evolution of cosmological per- 
turbations on linear scales. As discussed in [16], there is 
not a unique choice of these functions, and we have pre- 
sented results for two alternative choices: the pair used 
in [1], with (/i, 7) encoding deviations in the Poisson and 
anisotropy equations, respectively, as well as the pair 
(/i, E), with E in Eq. (6) directly related to the WL poten- 
tial. The main benefit of using these functions is that they 
allow for a model-independent test of the growth dynam- 
ics on cosmological scales even though they do not neces- 
sarily have a simple form in specific models of MG. Quite 
generally, in fact, they are defined through solutions of 
the equations of motion and depend on the choice of the 
initial conditions; still one can store observable informa- 
tion in these functions in a model independent way, and 
simply use care when translating the findings on fi and 
7 into results on the parameters of specific models [36]. 
In Sec. II B we have given a detailed description of how 
to perform the two-dimensional PCA of these functions 
with the aim of offering a useful technical reference for 
anyone who wishes to apply PCA to modified growth. 
In the same spirit, we have reviewed the observables and 
surveys used for our analysis in Sec. III. 

The bulk of the results is presented in Sec. IV, where 
we have analyzed the principal components (eigenmodes) 
of the MG functions for the combination of WL survey, 
CMB and SNe experiments detailed in Sec. IIIC. As al- 
ready noted in [1], the number of well-constrained eigen- 
modes gives a forecast of how many degrees of freedom 
describing deviations from GR will be constrained, and 
is particularly informative when comparing the outcome 
for different surveys. The shapes of the eigensurfaces in- 
dicate the regions of scale and redshift where the surveys 
under consideration will be most sensitive to departures 
from GR. 

We have given a detailed presentation of the eigen- 
modes and eigenvalues of the functions /i, 7, E as well as 
of the combination (/i, 7), comparing them and interpret- 
ing the differences. We have studied thoroughly the de- 
generacy between the MG functions and other cosmolog- 
ical parameters by progressively varying and marginaliz- 
ing over the different parameters. At every step we have 
interpreted and explained the effects of the marginaliza- 
tions. Of particular interest is the analysis of the degener- 
acy with the equation of state. We have found that after 
marginalizing over the MG functions, the high-z infor- 
mation on w(z) is erased and its best constrained eigen- 
modes are squeezed towards low redshift; however, as we 
show in Sec. VI, the effects of the marginalization are not 
dramatic and future surveys will have the ability to con- 
strain both w and modified growth. From the comparison 
of the results for LSST and DES, in Sec. IV E, we notice 
that LSST will have overall a higher sensitivity to modi- 
fied growth and will be more sensitive to time-dependent 
features. In Sec. V, we studied the effects of WL sys- 
tematics for LSST, and found that the degradation of 
constraints on MG is not significant, at least of the sys- 
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tematics models we have considered, and especially after 
one marginalizes over an arbitrary w(z). 

Finally, we have shown the utility of the PC A approach 
as a data compression stage. One can store the informa- 
tion contained in observables in terms of the MG pixels, 
or the eigenmodes of two functions. One can later use 
this information to project on to constraints on the pa- 
rameters of specific models. For example, in Sec. VII, we 
projected the errors on the MG functions to forecast the 
error in Bq. We have shown that only a fraction of the 
total number of eigenmodes is needed to obtain this con- 
straint. 

The degeneracy between [i and 7 or E can be further 
broken by adding information from redshift space dis- 
tortions (RSD) measurements [15]. It will be interesting 
to extend our analysis to include RSD for surveys that 
simultaneously measure lensing and peculiar velocities, 
such as DESpec and Euclid. 

Another direction for future work is to revisit the as- 
sumed scale-independence of the galaxy bias. In [57], 
it was shown that a scale-dependent growth necessar- 
ily implies a scale-dependent bias on linear scales. In the 
present analysis, this was effectively encoded in our func- 
tion /i, but in a future study it may be interesting to 
include the scale-dependent bias explicitly. In addition 
to scale-dependent bias, scale-dependent initial condi- 
tions are also likely to be degenerate with the MG modes 
(and probably worse, as they impact lensing as well.) Ar- 
guably, this is a more likely degree of freedom than the 
MG scale dependence. 

Our technique represents a model-independent way of 
analyzing the power of cosmological surveys to constrain 
modified growth. In addition to forecasts, it can be ap- 
plied to current data. For instance, given an array of 
experiments, one can use a Fisher forecast to first find 
the eigenmodes, then fit the amplitudes of these modes 
to real data. Because these modes are expected to be 
(nearly) orthogonal to each other, it does not matter if 
one fits them one by one or simultaneously. If any of 
them is found to deviate from zero significantly, it would 
constitute a smoking gun for modified gravity. Alterna- 
tively, if Fisher forecasted eigenmodes are found to be 
non-negligibly correlated, it would indicate that the fidu- 
cial model assumed in their derivation was wrong and 
that a modification is needed. One may also attempt 
reconstructing /i(fc, z) and 7(fc, z) from data, using the 
correlated prior technique introduced in [55, 58] for re- 
constructions of w(z). 
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Appendix A: The model for WL and GC systematics 

We follow [43, 44] and consider three sources of sys- 
tematics for future tomographic imaging surveys: photo- 
z errors, as well as additive and multiplicative errors due 
to the uncertainty of the point spread function (PSF) 
measurements. 



a. Photo-z Errors 

The redshift errors may stem from three sources: the 
distortion of the total galaxy distribution, z-bias and z- 
scatter. Suppose the distribution of the total galaxy is 



N(z) oc z 2 ex.p(-(z/z ) 2 ) 



(Al) 



where z denotes the true, spectroscopic redshift 
(throughout we use the over-bars to denote the precisely 
measured quantities to any precision), and zo is the me- 
dian redshift for a given survey. Then the galaxy distri- 
bution of the zth bin is, 
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(A2) 



where cr(z) denotes the photo-z error at redshift z, and 
erfc is the complimentary error function. 

Distortion of the total galaxy distribution: Suppose 
the observers measure the redshift z using some photo- 
metric method, and Az = z — z denotes the error. This 
error, in general, might induce a distortion of the overall 
distribution of the galaxies. To quantify this effect, we 
follow [44] to expand Az using the smooth Chebyshev 
polynomials (T n (x) = cos(narccos x)): 



N c hh - ley 
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As argued in [44], choosing iV c hb = 30 is large enough to 
yield convergent result, and we follow this setting. If the 
expansion coefficients g^s are much smaller than unity, 
the biased galaxy distribution can be estimated as 
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The galaxy distribution of the ith bin is, 
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(A5) 



Thus if the overall distribution is biased by g, all the red- 
shift bins are biased accordingly, as shown in Fig 22 (We 
show the galaxy distributions of LSST for an example). 
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FIG. 22: Redshift error - distortion: distortion of the total 
galaxy distribution due to uncertainties in redshift measure- 
ment, plotted here for LSST. g^s are expansion coefficients of 
Az in terms of Chebyshev polynomials (A3). 

In our calculations we account for this effect by 
marginalising over 30 Chebyshev coefficients. 

Redshift-bin Centroids Uncertainty: To be general, we 
also consider the possible degradation if there exists some 
uncertainty in measuring the centroids of the redshift 
bins, i.e. the so-called non- vanishing z-bias bi=Zi — Zi ^ 0. 
The z-bias basically 'shifts' the centroids of the bins. See 
this effect illustrated in Fig 23. In the calculation, we as- 
sign one z-bias parameter for each bin, then marginalize 
over them. 

Note that the centroid shifts do not capture the 
catastrophic errors where a smaller fraction of redshifts 
are completely mis-estimated and reside in a separate 
island in the z — z plane [44, 54]. 
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FIG. 23: Redshift error - shift: uncertainty in measuring the 
centroids of the redshift bins (non- vanishing z-bias), quanti- 
fied by bi coefficients at redshift bin i. 

z-scatter: We assume that <j(z) = ctq(1 + z), and we 



choose a = 0.03(0.05) for fiducial model for LSST 
(DES). But if do is not perfectly measured, i.e., the z- 
scatter is not zero, A = <Jo — (To ^ 0, there will be further 
degradation. In the calculation, we assign one z-scatter 
parameter for each bin, then marginalize over. As shown 
in [44], this effect is sub-dominant. 
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FIG. 24: Redshift error - scatter: uncertainties in the photo-z 
error at redshift z modeled as a(z) = cfq(1 + z), where Aj is 
do — (To at redshift bin i. 



b. Additive errors 

Additive errors are present for both galaxy counts and 
lensing shear measurements, and they are generated, for 
example, by the anisotropy of the point spread function 
(PSF). Following [44] and [43], we parametrize the addi- 
tive errors as 

(C? Y h=5xY P A?Aj(^y (A6) 

and choose p = 1,7 = 0. The fiducial values of the A's 
are, (A g ) 2 = lO" 8 , (A 1 ) 2 = 10~ 9 [43]. 



c. Multiplicative errors 

Multiplicative errors in measuring shear can be intro- 
duced by various sources. For example, a circular PSF of 
finite size is convolved with the true image of the galaxy 
to produce the observed image, and in the process it pro- 
duces a multiplicative error [44]. 

(C]) ij = (C]) ij [l + f i + f j ] (A7) 
We choose / = as fiducial. 
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